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A simple model of Laplacian growth is considered, in which the growth takes place 
only at the tips of long, thin fingers. In a recent paper, Carleson and Makarov used 
the deterministic Loewner equation to describe the evolution of such a system. We 
extend their approach to a channel geometry and show that the presence of the side 
walls has a significant influence on the evolution of the fingers and the dynamics of 
the screening process, in which longer fingers suppress the growth of the shorter ones. 
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I. INTRODUCTION 

A variety of natural growth processes, including electrodeposition, viscous fingering, so- 
lidification, dielectric breakdown or even growth of bacterial colonies can be modeled in 
terms of Laplacian growth. In this model a plane interface between two phases moves with 
velocity driven by a scalar field u{j,t) satisfying the Laplace equation 

V\(r,t) = 0, (1) 

with the boundary condition u{y, t) = at the phase interface. The normal velocity of the 
growing phase is proportional to the field gradient at the interface (or to some power {rj) of 
the gradient) 

|VM(r,t)r. (2) 

The harmonic field u can represent e.g. temperature, pressure, or concentration, depending 
on the problem studied. An important property of Laplacian growth processes is the Mullins- 
Sekerka instability of the advancing interface. The field gradient over a small bump is larger 
than that over the plane interface, thus, for ?7 > 0, the bump grows faster than adjacent 
areas of the interface and develops into a finger. Contrastingly, for < 0, the bumps are 
flattened and the growth is stable. 
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FIG. 1: Fingering in the combustion experiments fS] and the theoretical model (right) 



Initial phases of the evolution of a plane interface are well understood in terms of linear 
stability analysis, which yields the wavelength of the most unstable perturbation 
However, the later stages of the evolution are no longer linear and hard to tackle analytically. 
Here we consider a simplified model of the developed nonlinear state. It is assumed that a 
number of finger-like protrusions were already formed (as a result of the initial instability of 
the front) and the further growth of those fingers is assumed to take place only at their tips 
(see Fig.[l]), with velocities proportional to the field gradient. This dynamics is deterministic: 
once the initial geometry is given, the state of the system at any later time is uniquely 
determined. Thus, such a model can be used when, except for the initial instability leading 
to the finger formation, the role of a noise in the evolution of the system can be neglected. 
Additionally, we neglect another noise-driven phenomenon: the tip-splitting effect when 
single finger bifurcates into two or more daughter branches. 

The above model of finger growth was formulated mathematically by Carleson and 
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Makarov [4| (see also Selander's thesis [5|) and called by them „the geodesic Laplacian 
path model". Independently, a similar idea was considered by Hastings in j^. 

In several experimental and numerical studies the patterns relevant to the above- 
introduced model have been observed. The examples include dendritic growth in some 
of the electrochemical deposition experiments , channeling in dissolvin g ro cks [8| , side- 
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branches growth in crystallization [9|, [l^, or fiber and microtubule growth 
the most beautiful experiments on the fingered growth are the combustion studies by Zik 
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and Moses [3|, Il3|. In those experiments (some of the results of which are reproduced in 
Figs. [H and [2]) a solid fuel is burnt in a Hele-Shaw cell, i.e. in the narrow gap between two 
parallel plates. Near the fiame extinction, as the fiux of the oxygen supplied to the system is 
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FIG. 2: Competition in the fingered growth in combustion experiments [3!|. Initially equal sized 
fingers (bottom) evolve towards state where every other finger stops growing (top). 



FIG. 3: Competitive dynamics of the channels in the dissolving rock fracture. The figures present 
the dissolution patterns at two different time points. Similarly to Fig.[2l longer channels grow faster 
and suppress the growth of the shorter ones ^. 

being lowered, the initial instability of the combustion front develops into the sparse fingers 
(cf . Fig. [l|) , which appear to evolve in a regular or near-regular way. 

Growth of the fingers demonstrates analogous instabilities as the plane front, since the 
gradients around the tips of the longer fingers are larger than those around shorter fingers. 
This leads to so-called "shadowing effect" - the longer fingers grow faster and suppress the 
growth of the shorter ones in their neighborhood, which in some cases gives rise to a scale- 



invariant distribution of finger lengths in the long-time limit [8|, uA, ll5|. Figs. [2] and [3] 



present examples of such a competitive growth in the channeling processes in dissolving 
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rock and in the combustion experiments described above. 

Most of the experiments on nonequihbrium growth mentioned above were performed in 
a quasi-2d geometry. A convenient way of solving the Laplace equation in two dimensions is 
to use a conformal mapping which transforms a domain under consideration to some simpler 

n 

region where the solution may easily be found. A remarkable idea, due to Loewner [16|, is to 
trace the evolution of the mapping instead of the evolution of the boundary itself. It turns 
out, namely, that the evolution of the map may be in many cases described by the first order 
ordinary differential equation (Loewner equation), which represents a considerable simpli- 
fication in comparison to the partial differential equation describing boundary evolution. 
Loewner evolutions are intensely studied in the theory of univalent functions (for general 
references see the monographs l3, [3| and for a recent physical introduction see Q, S|)- 
The subject has recently attracted a lot of attention in the statistical physics community in 
context of Stochastic Loewner Evolution (SLE), which has proved to be an important tool 
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in the study of two-dimensional critical systems 

The exact form of Loewner equation depends on the shape of the domain in which the 
growth takes place. Usually, it is either the complex half-plane (where the initial phase 
boundary corresponds to the real axis) or radial geometry (where the initial boundary is the 
unit circle in the complex plane). However, many experiments on the unstable growth are 
conducted in a channel geometry, between two reflecting walls. In this paper we show how 
to extend the Carleson and Makarov model to such a geometry. Using conformal mapping 
formalism we derive Loewner equation for that case, which allows us to flnd the dynamics 
of the fingers and analyze the shadowing process. As it turns out the influence of the walls 
is often crucial for the dynamics of the fingers. 

II. THE MODEL 

With this introduction, let us formulate the model of finger growth to be considered. The 
growth takes place at the tips of a finite number, n, of infinitely thin fingers Ti(t) (disjoint 
Jordan arcs) 

r,(t)cW, i = l,...n, (3) 
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where W is the domain in which the growth takes place and where the Laplace equation 
needs to be solved. The fingers extend from the boundary of W towards its interior in such a 
way that Ti{t) C Ti(t') for f > t. Additionally, on both the phase interface (along boundary 
of W ) and along the fingers the condition u(r, t) = is imposed. 

Since the finger is assumed to be infinitely thin, there is a singularity in a field gradient 
at its tip. Namely, at a small distance r from the tip of zth finger, the gradient takes the 
form 

VM(r, t) = ^ (cos(e/2)e, + sin(^^/2)ee) , (4) 

where the coefficients Ci{t) depend on lengths and shapes of all the fingers. In the above, 
the origin of coordinates is located at the tip of the finger and the polar axis is directed 



along it. Following Derrida and Hakim 



25| . we introduce a small circle of radius tq around 



the tip and define the finger growth rate as the integral of field gradient over the circle 

v,{t) = in- VM(r, t) ds = 2^ Ci{t). (5) 



Note that if the field u{j, t) describes the concentration, then the above integral corresponds 
to the total particle fiux through the circle. The parameter tq should be of the order of the 
finger width; its exact value does not infiuence the dynamics as long as we assume it to be 
the same for each finger. In such a case, the factor may be absorbed into the definition 
of time, and we take Viit) equal to Cj(t) (or to Ci(t)'^ in r] growth). 
To solve the Laplace equation we construct a time dependent map gt: 

gt:W\{Ti{t)U...UTn{t))^W (6) 

together with its inverse, ft: 

w^w\(ri(t)u...ur„(t)), (7) 

ftogt = gto ft = id , t>0, (8) 

as illustrated in Fig. [H The function ft can be extended to a continuous function of the 
boundary, which is two-to-one along the fingers, except of the tips where it is one-to-one. 
Thus the tips of the fingers (denoted by 7j(t)) may be added to the domain of the function 
gt with the corresponding images, aj(t): 

fli = 9t{7i)- (9) 
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FIG. 4: An example configuration of three fingers (ri,r2,r3) extending from the real axis in the 
physical plane (z-plane). The mapping gt maps the exterior of the fingers onto the empty half-plane 
(w-plane). The images of the tips 7^ are located on the real line at the points x = Oj. The gradient 
lines of the Laplacian field (dot-dashed) in the z-plane are mapped onto the vertical lines in the to 
plane. In a given moment of time, the fingers grow along the gradient lines the images of which 
pass through the points a,. 

It may be shov^n (see Appendix [A] and also Refs. ay , l25|) that the growth rate of a 
finger ([5]) can be expressed in terms of the map ft as 

^;,(t)~|/;'(a,(t))|-''/^ (10) 

As mentioned in the Introduction, the above model of finger growth was proposed in 



Refs. 0, 0] and independently in 



6|. This model may be also looked upon as a deter- 



26l . I27I I needle growth model - a simplified, 



ministic generalization of Meakin and Rossi 
non-branching version of Witten and Sander's diffusion-limited aggregation (DLA) [28|. In 
Meakin and Rossi model the growth is allowed to occur only in the direction away from the 
substrate, which results in a forest of parallel needles with a broad distribution of heights. 
Gates 29j obtained average density of the aggregate in such a system in mean field approx- 



301 calculated fractal dimensions of needle structures for different 



imation whereas Evertsz 
values of the exponent rj. 

Yet another modification of DLA was studied in so-called "polymer-growth model" 
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where, as before, the particles were allowed to attach to the tip of the growing chain 
only, but the condition that the growth must occur in the direction away from the substrate 
was relaxed. This time, the process resulted in a set of fiber-like chains of different lengths 
and shapes. 
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The above described, DLA-based models are stochastic in nature. On the other hand, the 



deterministic versions o: 
authors 
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le Meakin and Rossi needle growth model were studied by several 
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37, 
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39l |. mostly by conformal mappin g^ te chniques, which 



were first applied to DLA-like aggregates by Shraiman and Bensimon Ball [4ll|, and 
Szep and Lugosi [42|. The majority of authors considered the problem in a radial geom etry, 
with a set of straight needles growing radially from the origin, whereas in a set 

of parallel needles with alternating lengths was analyzed in periodic boundary conditions. 
An important difference between the above quoted needle models and the model considered 
here is that we do not constrain the fingers to follow the straight line, instead they can 
bend in the direction of the field gradient at the tip. Thus, whereas the needle models 
are well-suited to describe a strongly anisotropic growth of rigid structures such as the 
sidebranching dendrites in solidification, they are less suited to describe the phenomena in 
which the growing structures may become defiected by the field, as it is the case in the 
above-described combustion experiments, channel formation in the dissolving rock or in a 
number of fibril-growth processes. 

In the next section, we sketch the derivation of the Loewner equation for the fingered 
growth in the half-plane. In the context of the Carleson-Makarov model, this problem was 
analyzed by Selander in j5|. We repeat this derivation here, because it is the simplest case 
which illustrates the method involved, which will be then used to tackle the channel geometry 
case. 



III. FINGERED GROWTH IN THE HALF-PLANE 

In this section, following we present the description of the fingered growth in the 
complex half-plane in terms of the conformal mappings. The solution of the Laplace equation 
in the empty half-plane, satisfying the boundary conditions: u{r) = on the real axis and 
1^ — > 1 ( constant fiux at infinity), is given simply by 

y ?/— >oo 

u{x,y)=y. (11) 

The map gt from the exterior of the fingers to the empty half-plane will be uniquely deter- 
mined by the so-called hydrodynamic normalization at infinity 

lim gt(z) -z = 0, (12) 



8 



which ensures that the flux at infinity is unaffected by boundary movements. One of the 
possible ways of finding the map gt is to construct it as a composition of slit mappings, with 
the slit length going to 0. In the case of the upper half-plane C+ = {z E C : Re{z) > 0}. 
the slit mapping 

: C+\(a + z[0,/i]) ^ C+, h>0 (13) 

is of the form 

(j){z) = ^{z-ay + h^ + a. (14) 

In Eq. ffT4l) the square root is specified by demanding that (pi^z) — > z at infinity, thus the 
difference (0(x) — a) is negative on the real axis for x < a and positive for x > a. As stated 
above, there is a branch cut along {a + i [0, h]); the left and right side of this segment are 
mapped on the interval [—h, 0] and [0, h] respectively. 

In the simplest case of the single finger (n = 1) the idea of the construction of the mapping 
by the composition of successive slit maps is shown in Figure [5l The map gt maps the arc 

|r(f + r) — > ^ 

ait) 



a(t + T) 

FIG. 5: Illustration of the composition of conformal maps described in the text. The mapping gt+r 
is obtained as the composition of gt and the elementary slit mapping (fT4l) . 

T(t) into the real axis R, while gt+r maps the longer arc Tit + r) into R. The image of the 
tip, a, is also in general time-dependent, a = a{t). For small r, the mapping gt+r may be 
obtained by the composition 

gt+r = (l>ogt + 0{T^), (15) 

where (j){z) is the above defined slit mapping (fT4ll . Since the growth rate of the finger is not 
constant, the slit length, h, will be a function of time. It is convenient to choose h in the 
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form 

h = ^/2Td{t), (16) 

where d{t) is the so-called growth factor [5|. The square root dependence of the slit length 
on the timestep ensures that the increment of length of the finger will be linear in r. The 
slit mapping (f){z) is thus given by 



= ^/{z - a(t))2 + 2rc/(t) + a(t). (17) 



In the composition (iTSl l only the terms linear in r are needed, hence 

g,^, = (Pog, = g, + -^^^- + 0{T^), (18) 
9t - a{t) 

^^±^ = ^^ + 0(r). (19) 
r gt- a{t) 

In the limit r — we obtain the Loewner equation for one finger in the half-plane 

* = (20) 

gt - a{t) 

with the initial condition 

go{z) = z (21) 



corresponding to the empty space with no fingers. Note that the pole of the equation ( l20ll 
is located at the image of the tip, a{t) = gti^l). To relate the growth factor d{t) to the finger 
velocity let us notice that for a finite r, the composition (p o gt increases the length of the 
finger approximately by Td{t) |/"(a(t))|, thus the growth velocity is given by 

v{t) = d{t)\f:{am- (22) 

Hence, using f fTOl ) we get 

d{t) = \f:{a{t))r'"-"- (23) 

The final element of our description is the position of the pole a as a function of time. This 
function may be found from the condition that the finger grows along the direction of the 
gradient near the tip. Due to the singularity at the tip, it is more convenient to work in the 
C(j-plane. Namely, the gradient lines in the physical plane (z-plane) are mapped by gt onto 
the vertical lines in the a;-plane, as illustrated in Fig. [H The counterimage of the veritcal 
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line uj = a + is, s > will thus define the growth direction in the z plane. In the case of 
a single finger, from symmetry, the gradient line in the z plane is also a vertical line - the 
one passing through the tip of the finger. Thus in this case the position of the pole must be 
constant 

a{t) = const. = a(0). (24) 

The Loewner equation can be generalized to the n-finger case by analyzing the composi- 
tion of n slit mappings, (pi, one for each finger. Analogously to (fTTl) . the slit mappings are 
then given by 

(p^{z) = ^{z - ai{t)f + 2rdi{t) + ai{t), (25) 

where aj(t) and diif) are the image of the tip and the growth factor of the ith finger, 
respectively. This leads to the Loewner equation of the form 0, 5| 

9i = y^^,- (26) 

This time, however, to force the fingers to grow along gradient lines, the functions aj(t) need 
to change in time, since the tip images Oi = g{^i) are moved by slit mappings 0j with j ^ i. 
This leads to the following condition for the motion of the poles 

Let us note at that point that the above idea of generating the growing aggregate by 
iterated conformal maps was applied also to the original DLA problem in a seminal paper 



by Hastings and Levitov [43||. In fact, the "strike'-mapping proposed by them in [43!| is a 
counterpart of the slit mapping ( fT4l l in the radial geometry. Analogous constructions for 
the DLA in a channel and cylindrical geometry were proposed in Refs. Q, [4^. However, 
unlike the model considered here, those models were stochastic in nature, and - as usual 
in DLA - generated noise-driven, branched, fractal structures. A deterministic version of 
Hastings and Levitov construction was proposed by Hastings in [^. Although he did not use 
the formalism of the Loewner equation, his model is essentially analogous to that presented 
above. 

Going back to the system of equations (I26II27I) . let us now look at the solutions in a 
few simplest cases. As mentioned, the single finger case is rather straightforward: the 
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finger grows vertically along the line x = a. A more interesting case is that of two fingers. 
Here the results depend in a significant way on the value of the exponent rj used. A related 
problem in a slightly different geometry (growth in W = C\R+) was considered by Carleson 

n 

and Makarov in [J]. Although their results are not directly applicable here, it is relatively 
straightforward to repeat their derivations for the growth in the half-plane. Namely, it turns 
out that there are three regimes in the behavior of the fingers, depending on the value of 
the exponent t]. For — oo < r] < rjc ^ 0.43, the two fingers, irrespectively of their initial 
positions, grow symmetrically, as illustrated in Fig. [6l Next, for rjc < r] < r]'^ = 2/3, the 
symmetric solution becomes unstable; one of the fingers starts to grow faster and screens 
the other. However, the screening is only partial and the ratio of finger velocities Vi/v2 goes 
to a positive constant (different from 1) when t — oo. Finally, for r] > t]'^ there is a stronger 
screening and the ratio of the velocity of the slower finger to that of the faster one goes to 
zero asymptotically. 

1.5 
1. 

y 

0.5 



-0.5 -0.3 -0.1 0.1 0.3 0.5 -0.5 -0.3 -0.1 0.1 0.3 0.5 

X X 

FIG. 6: Two fingers growing in the half-plane for r] = —2 (left panel) and r] = A (right panel). The 
initial positions of the fingers are ai(0) = —02(0) = 0.1. 

The shape of the fingers in the first regime may be obtained analytically. To this end, it 
is most convenient to consider the case off] = —2 which, according to Eq. f l23ll corresponds 
to the evolution with constant growth factors, di = d. For further analysis, we choose the 
coordinates in such a way that ai(0) = —02(0) = a^. Eq. (j27|) may be then readily integrated 
to yield 

ai,2(t) = ±Jal + dot. (28) 
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Inserting the above into the Loewner equation (l26ll . we can find gt and then obtain the 
positions of the tips 7i(^) cis a function of time. They are given imphcitly by 

± 16{al + dotf/^ = 7i,2(t)(7i,2(t)' - 5aiy, (29) 

where we take only the roots which lie in the upper half-plane and satisfty the initial con- 
dition 7j(0) = aj(0). An example trajectory with ao = 0.1 is shown in Fig. [H We see that 
the fingers repel each other, which is due to the term l/{aj — ai) in the evolution equation 
of the poles ( 1271 ). Asymptotically, as t — * oo, the relation f l29l ) takes the form 

± lQ{dotf'^ = (7.(t))^ (30) 

from which it may be concluded that the fingers tend to the straight lines arg(z) = and 
arg(z) = ^ with the angle 7r/5 between them. Although the above result was obtained for 
r] = —2 case, as mentioned above, such a symmetric solution remains stable up to rjc. Thus, in 
this range of i] the fingers follow the same trajectories as in Fig.[6l but with different velocities 
than those in d = const, case. In fact, whenever the fingers grow symmetrically, with 
equal growth factors, di(t) = d2(t), one may always rescale the time coordinate by defining 
t' = di{t")dt" and thus reduce the problem to the constant growth factor (di = ^2 = 1) 
case. 

The right panel of Fig. [6] shows the shape of the finger in the strongly unstable case, 
?7 = 4. In this case the solution was obtained numerically, as described in Appendix [Bl As 
it is observed, initially the tips follow the same trajectory as in the stable case, but then 
due to the numerical noise in the computations, the instability sets in, one of the fingers 
outgrows the other and then continues along the vertical direction. The other finger slows 
down and the ratio of its velocity to that of the winning finger goes to zero. 

The fact that there is a range of positive values of r] for which the symmetric solution is 
stable seems surprising at first sight. However, in the unbounded domain the fingers have 
the possibility to escape from each other to the regions where the infiuence of another finger 
is smaller. Such a behavior of the fingers may be hardly observed in a bounded system, 
because of the wall effect which makes it impossible for the fingers to escape from each 
other. In the next section, we show how to take into account the presence of the side walls 
in the system and how their presence affects the dynamics of the fingers. 
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FIG. 7: The geometry of the channel together with the boundary conditions 
IV. FINGERED GROWTH IN THE CHANNEL GEOMETRY 

In this section we consider the growth of the fingers in the channel geometry, i.e. in the 
domain 

V = {z = x + iye€:y>Q, a;G]-l,l[}, (31) 

with the boundary conditions u = at a bottom wall [—1, 1] and along the fingers, and g = 
on the sides (which corresponds to reflecting boundary conditions at the impenetrable side 
walls) and |^ = 1 at infinity (see Fig. [71). For an empty channel with the above boundary 
conditions, the solution of Laplace equation is again given by 

u{x,y)=y. (32) 

in full analogy to the half-plane case (fTTl) . 
Our goal is to construct the mapping gt 

(7t:P\(ri(t)U...ur„(t))-.P (33) 

together with its inverse ft- The boundary condition at the bottom wall is different from 
that at the sides, thus the points —1 and 1 must remain fixed under the mapping ft, or, in 
terms of Qt 



lim gt{z) + 1 = \im gt{z) -1 = 0. 

z—>—l z— >1 



(34) 
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Additionally, we require that gt keeps the point at infinity fixed, i.e. 

lm{gt{z)) — > oo for lm.{z) — > oo (35) 

Those conditions define the conformal map gt uniquely. 

First let us consider the case of the single finger. The derivation of Loewner equation 
in this case will be analogous to that presented above for the half-plane. However, the slit 
mapping is now different. The idea of construction of the mapping is shown in Fig. [3 As 
before, we begin with a short slit of length h = ^Ird. The first function '1' transforms P 
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FIG. 8: Schematic picture of the mappings used in construction of elementary slit mapping (j) in 
the channel geometry 

into the half-plane keeping points A = —1, E = 1 fixed and is given by sin (f-z). When the 
slit length is small, h « 1, its length after the transformation T can be calculated from 
the Taylor expansion up to linear terms in r 

/ \ / TT^ \ / TT \ ITT /TV \ 

sin i^-{a + iV2T^) j ^ i 1 + —Td\ sin i^-aj + Y^^^^ ■ (36) 

Denoting the new length of the slit by V^rd' we get 

2 

d' = d^cos^ (|a) +0(r2). (37) 
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Function '2' is the slit mapping in the half-plane given by (lT4ll . The composition of '1' and 
'2' reads then 



^{z) = ^ (sin {^z^ -(3 sin (|a)) + 2rc?', (38) 

where /? = (l + ^Td^ . The above slit mapping moves the points A and E. To shift them 
back to -1 and 1, we use additional map, '3', which is a linear function with real parameters. 
The map '4' is just the inverse of '1' , i.e. | arcsin(2;). The final form of (p is thus given by 



(z) = — arcsin 

TT 



Keeping only the terms linear in r 



2v.(z)-(y(l) + V9(-l)) 



(39) 



z sm [^z) — sm [^a) 

we obtain Loewner equation of the form 

,/ N TT cos i^gt) , . , 

9t = ■ X % ,,,v (41) 

2 sm (^fi-tj - sm (2a(t)j 

Note that due to the presence of the side walls the slit mapping is no longer symmetric, in 
the sense that the images of the points B and D at the base of the finger are asymmetric 
with respect to the image of C (cf. Fig. [8|). This means that, in contrast to the half-plane 
case, the pole a = g{j) will be shifted by the mapping. This shift may be obtained from the 
Loewner equation in the limit g{z) a. However, the equation is singular at that point, 
thus we take the symmetric limit from both sides towards a singularity: 

, , , Wia - e) + Wia + e) 

€— >o z 

where W{g) = cos{^g)/ (sin(|5f) — sin(|a)). Explicit evaluation of the limit yields 

a(t) = -^d{t) tan (^a(t)) . (43) 

The equation for the n— finger case may be obtained, as before, by the composition of n slit 
mappings, one for each finger, which leads to 
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Finally, to derive the condition for the motion of the poles in the n— finger case, we 
need to add self terms of the form (l43l) and the interaction terms, which may be obtained 
from ( !44l ) by taking gt = aj with i j. This leads to 

cos (f flj) 



i = l 



sin (f aj) — sin (f aj) 



(45) 



The presence of the self term, ( l43ll . attracts the pole to a = 0, which is a stable fixed point 
of Eq. (l43l) and causes the finger to bend in the direction of the centerline of the channel. 
It is most clearly seen in the single finger case. For the constant growth factor, d(t) = do 
(which corresponds to r] = —2) equation (l43l) may be solved explicitly to yield 



"K 

sm I -at 



TT 

sm I -ao 



(46) 



Thus, sin(|at) goes to zero exponentially. A corresponding trajectory of the finger may be 
expressed implicitly through elliptic integrals and is shown in Fig. [9l As it is observed, the 
finger starts at z = a(0), initially grows perpendicularly to the bottom wall, but very soon 
the infiuence of the walls becomes important and the finger is attracted to the symmetric 
position in the center of the channel. As explained previously, for other values of rj the shape 
of the finger is the same as that presented above, only its velocity changes. Naturally, if 




FIG. 9: The growth of the single finger in the channel with a(0) = —0.5 



the initial position of the finger is already in the middle of the channel, a(0) = 0, it would 
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simply continue growing along the centerline. In that case the map ft is given by 

Mz) = ^ arcsin (^^sin^ cosh^ - sinh^ (^^(^))) ' (47) 

where H{t) is the height of the finger at a given moment of time. In particular, 

|/;'(0)| = ^coth(|/J(t)), (48) 
which gives the velocity of the tip as 

v{H) = [^coth {^H)y'^\ (49) 
The asymptotic velocity of the finger is thus 



= \im^v{H) = i^-) , (50) 
and the asymptotic growth factor is (cf. Eq. (i23ll ) 

das = (-) . (51) 

It is instructive to compare this result with that for a single finger growing in the half plane. 
In the latter case 



and the velocity obeys 

v{H) = H^'\ (52) 

hence, for 77 > 0, it is growing indefinitely as the finger increases its height. This is to 
be expected since in the unbounded the finger is getting further away from the 

absorbing wall, it intercepts an increasingly larger fiux. Contrastingly, in the case of the 
channel, the total flux through its cross-section is finite, and the growth rate of the finger 
stabilizes as soon as its height becomes large in comparison to the channel width. Note that 
the result (l52l) for the half-plane may also be recovered by expanding the formula ([49l) in 
H <ti I, since | coth (f i?) = H + . . . . Thus, at the beginning of the evolution, the finger 
behaves as if it were in an unbounded domain; soon, however, the presence of the side walls 
becomes a determining factor in its dynamics. 
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V. FROM THE CHANNEL TO THE CYLINDER 

The above formalism may also be used to find the evolution of the fingers in a channel 
with periodic boundary conditions: 

u{x + 2,y) = u{x,y), (53) 

topologically equivalent to the surface of a semi-infinite cylinder. Namely, consider a sin- 
gle finger growing in the periodic channel with the initial position of the pole a(0) = oq. 
Translating the origin of the coordinate system to we obtain a single finger growing ver- 
tically in the middle of the channel on the side walls of which both periodic and refiecting 
boundary conditions are satisfied simultaneously. The growth of this finger is described by 
the Loewner equation (1411) with a = 0. Transforming back to the original coordinates one 
obtains 

gt = cot (^^{gt - a)) , (54) 

which is the Loewner equation for a single finger growth in the cylinder. This equation has 
already been derived j^, 20| in the context of growth processes in radial geometry (which may 
be then mapped onto cylindrical by the map of the form z log 2;). However, yet another, 
and perhaps easier, derivation of (l54ll may be given, starting form the Loewner equation for 
the half-plane (|20|) and summing over the periodic images of the pole: aj = a + 2j, j G Z. 
This gives 

The generalization of the Loewner equation to the n— finger case proceeds along the same 
lines as before and yields 

gt = ^d~ cot (^^{gt - ai)^ , (56) 

i 

whereas the equation of motion of the poles reads 

% = {^^^^ ~ ■ ^^'^^ 

i 

This time the self term (l43ll is absent since the situation is again symmetric and the single- 
finger slit mapping does not affect the position of the corresponding pole. 
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FIG. 10: The channel A and its mirror reflection A' 



VI. FROM THE CYLINDER TO THE CHANNEL 



In the previous section we derived the dynamics of the fingers in the cylinder based on the 
Loewner equation for the channel. Here we will go in the opposite direction and derive the 
Loewner equation for the channel starting with Eq. (i54l ). As it turns out, such an approach 
provides us with the clear interpretation of the self term (l43ll . 

To start with, consider n fingers growing in a channel with refiecting boundaries (desig- 
nated by A in Fig. [101). Next, let us refiect the system with respect to one of its side walls 
and denote the image by A'. Thus, there are now 2n fingers: the original ones {1, . . . , n} in 
A and their images {n + 1, . . . , 2n} in A'. Due to the refiection symmetry 



i E {I. ..n}, 



(58) 



with coordinates chosen so that x = corresponds to the joint wall of the two channels, as 
in Fig.Hni Finally, we look for the solution of the Laplace equation in the joint system (AA') 
in periodic boundary conditions. Due to the symmetry of the problem, such a solution will 
automatically satisfy the refiecting boundary conditions in A. 

To be more formal, the Loewner equation in the AA' channel will be given by 

2n 

gt = ^^di cot (^^{gt - tti)^ , (59) 



i=l 



which is equivalent to (l56ll additionally rescaled by a factor of 2 to account for the fact that 
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the width of AA' is twice the width of A. Taking into account (l58ll . we get 



gt = - "^d, l^cot (^-{gt - ai)j + cot (^-{gt + 

i=l 



(60) 



which, by straightforward trigonometry, leads to 



9t 



TT 



sin {^gt) 



,=1 i^^gt) - cos (f fli) 



2 cos 



(61) 



Finally, we move the origin of coordinates to the center of A by 

z' = z+l, 

which leads us to the Loewner equation for the channel ( l44ll . 

The equation of motion of the poles may be transformed along the similar lines: 

2n 



(62) 



TT /TT 

— cot — a,- — a, 
4 V4^ ^ 



cot ( — (flj — fli) 



(63) 



The first term describes the interaction between the pole aj and its image a„+j and, us- 
ing (l58l ) and fl62l) . may be transformed to —\dj tan (fflj) , which is exactly the self term f l43ll 
derived before. The terms in the remaining sum, using ( l58ll . may be written in the form 
cos(|aj)/(sin(|aj) — sin(|aj)), in agreement with ( l45l) . 

The above derivation shows, in particular, that the stable, symmetric solution of two- 
finger growth in the cylinder may be reduced to the previously considered problem of a 
single finger growth in the channel. Indeed, Fig. [TT] shows two-finger solution in the cylinder 
the initial conditions ai,2(0) = ±0.25 (if the fingers are initially not symmetric with respect 
to zero they can be rendered so by an appropriate change of variables) . It is observed that 
the right-hand side finger has the same shape as that in Fig. [9l whereas the left-hand side 
one is its mirror refiection. 

Fig. [m illustrates an important property of stable solutions of n fingers growing in the 
cylinder. Namely, irrespectively of initial positions, the fingers end up in the symmetric 
configuration, with equal distances between each other. This is a direct consequence of the 
repulsion between the poles, as described by Eq. (I57D . Indeed, it is seen that the uniform 
arrangement of the poles constitutes a fixed point of f lSTl) . 

Because of the relation between the geometry of a refiecting channel and that of a cylinder, 
the above statement may also serve to prove that, in the stable case, both the fingers and 
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FIG. 11: Two fingers growing in the cylinder with 01,2(0) = ±0.25 and 77 = —2 (stable symmetric 
solution). 

the poles in the channel end up at the positions Xi = {2i — l)/n — 1 where n is the number 
of fingers and i = 1, . . . ,n. 

VII. THE COMPETITION BETWEEN FINGERS 

In this section we look more closely at the instabilities in the growth of two fingers both 
in the channel and in the cylinder, paying particular attention to the screening process and 
competition between the fingers. The results presented in this section, unlike those presented 
in Figs. [9] and [m are obtained by numerical calculation (as described in Appendix [Bll . 



-!Z(U-I.)la 




« = 



u ~ (IJl,)-'"' 



FIG. 12: The scaling of a Laplacian field u in rectangular and wedge-like fjords. 



As mentioned in Sec. IIIIl in the case where fingers grow in the half-plane, there exists a 



22 



critical value of the exponent (rjc) below which the symmetric solution is stable nad another 
threshold value, r/^, such that for rjc < r] < rj'^, there exists an asymmetric solution of the 
growth of two fingers with the velocity ratio of the slower finger to that of the faster one, 
^, asymptotically approaching the value different both from and 1. Here we argue that 
for the channel, either reflecting or periodic, rjc = Vc = 0? i-^- the growth is unstable for any 
positive 7], with ^ ^ 0. Such a fundamental difference between the two systems is connected 
with the geometries involved Consider the late stages of the screening process, with a 
shorter finger (of length li) situated in a deep fjord between two longer fingers (or a longer 
finger and its periodic image) of length I2. The value of the Laplacian field u near the tip of 
the short finger then obeys 

u ~ e-"('^-'^)/", k > h, (64) 

where a is the width of the fjord. As mentioned, in the half-plane case the fingers asymp- 
totically tend to straight lines growing radially from the origin, thus in that case the fjords 
will have a wedge-like shape (cf. Fig. [121). For such a fjord with an internal angle 9 one finds 

I \ --^1^ 



[fj ' ^2 > h. (65) 

Thus the screening in the channel geometry is exponential, whereas in the half-plane it 
follows a power-law. Consequently, in the latter case, we can choose the values of r] and 9 
in such a way that would guarantee that the ratio of the velocities of the fingers is equal to 
their length ratio 

= ^. (66) 

which corresponds to the (asymmetric) solution of the fingers dynamics which is stationary 
in the sense that the ratio h/h remains constant. 

Contrastingly, for the exponential screening, irrespectively of the values of a and r] (as 
long as it is positive) the velocity of the shorter finger always decreases to zero, thus also 
h/h 0: the screening is complete. 

Fig. [I3] presents an example of such a situation. Here, two fingers are evolving in the 
channel with either reflecting or periodic walls at 77 = 1. At the beginning, the flngers 
repel each other and behave very similarly to the stable, symmetric solution with 77 < 
(cf. Fig. [H]). However, as soon as the their height becomes comparable to the width of 
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FIG. 13: Two fingers growing in the cylinder (left panel) and in the reflecting channel (right panel) 
for r] = 1. The initial positions of the fingers are ai(0) = —0.2 and 02(0) = 0.4. Inset: the 
interaction of two fingers in the combustion experiments of Zik and Moses [3j|. 

the channel, the instability sets in and initially small differences in height of the fingers are 
rapidly amplified. As observed in Fig. [131 i^i the refiecting channel the longer (winning) 
finger is attracted to the centerline (as it is the case for the single finger solution), whereas 
in the cylinder it continues to grow vertically. 

The analysis of Fig. [13] together with the corresponding dependencies of the growth 
velocity on time presented in Fig. [TSl shows distinctly three different stages of finger growth. 
In the initial stage, the fingers repel each other, their velocities increase quickly but remain 
equal to each other; the competition between the fingers is not yet present. Then, as the 
distance between the fingers approaches 1, their velocities stabilize at around v ~ 0.6. In 
fact, the velocities of the fingers in this stage are close to the asymptotic velocity of a single 
finger in the channel of half the width of the original one, i.e., analogously to Eq. (jSOll 

/7r\ 

v= i^-J ~ 0.56, (67) 

which agrees with the value 0.6 quoted above. As seen in Fig.[T4]the unstable 77 = 1 solution 
diverges from the stable one {rj = —2) when the height of the finger reaches approximately 
1—2 units, i.e. becomes comparable with the channel width. This is the moment where the 
final, third stage of the evolution begins, accompanied by the sharp decrease increase of the 
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FIG. 14: Two fingers growing the reflecting channel at r/ = 1 (sohd hne) and rj = —2 (dashed). The 
initial positions of the fingers are the same as those in Fig [131 



velocity of the losing finger and the corresponding increase of the velocity of the winning 
one up to the asymptotic value given, according to ((501), by {7r/2)-^^^ ^ 0.798. 

Note that in the case of the cylinder it is impossible to predict beforehand which finger 
is expected to win the competition - the initial situation is fully symmetric and in principle 
the system could remain in the unstable symmetric state were it not for the numerical noise 
(no other sources of noise are present). Contrastingly, in the refiecting channel, the finger 
which starts closer to the centerline is destined the win, since, during its evolution it moves 
through regions of higher field gradient. This is also the reason why the instability sets in 
earlier in the case of the channel. 

Further insight into the competition process may be gained by analyzing the evolution 
of the poles, as presented in Fig. [161 It is seen that in the channel the pole of the winning 
finger is attracted towards the center, whereas the pole of the suppressed one moves to the 
side of the channel. Such a behavior of the poles may be explained based on the evolution 
equation f[45ll in which the growth factors di have been replaced by their asymptotic values: 
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FIG. 15: The velocities of the fingers as a function of time in the cyhnder (left) and in the channel 
(right), corresponding to the fingers of Fig. [131 
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FIG. 16: The positions of the poles as a function of time, in the cylinder (left) and in the channel 
(right), corresponding to the fingers of Fig. [131 

di = das (cf. Eq. (EH)) for the winning finger and d2 = for the losing one. This leads to 



71 fir 
ai = -^das tan (^-ai 



^ . (68) 

2 sm (g'^aj — sm (2*^1) 

The asymptotic fixed point of the first equation is ai = 0, in full analogy to the single finger 
solution in the channel. It is then inserted into the second equation to yield 

d2 = —-rdi cot ( — 0-2 ) (69) 



4 V2 

with the asymptotic fixed points at 02 = ±1. 

In the inset of Fig. [131 we reproduce a two-finger pattern from the combustion experiments 
by Zik and Moses. The resemblance is remarkable. Nevertheless, we do not claim that such 
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a simple model can accurately reproduce the complex dynamics of combustion in Hele-Shaw 
cell, we can only conclude that the fingers in those experiments indeed seem to follow the 
gradient lines of the Laplacian field, (which in this case represents the concentration of 
oxygen at a given point) and that there is a screening mechanism leading to the competition 
between the fingers. This does not mean necessarily that the velocity of the finger must be 
connected with the field gradient analogously to Eq. ((21). In fact, many of the results of the 
above presented model are independent of the actual relation between the velocity v and the 
field u. In particular, the shape of the single-finger solution and the symmetric two-finger 
solution are independent of v{u) relation, which affects in this case only the overall growth 
rate, not the finger shape. The same does not hold for asymmetric solutions, but even in 
that case some of the qualitative features of the dynamics, which give the winning finger in 
Fig. [13] its characteristic shape remain universal: at the beginning the fingers are attracted 
to the symmetric solutions, then the instability sets in and finally the winning finger is 
attracted to the centerline. 



VIII. SUMMARY 



In this paper we have studied a simple model of Laplacian growth, in which the finger- 



like protrusions grow only at their tips, 
previously in the literature, 



Q, HQ, 



n contrast to the needle-growth models, considered 



35|, 



37 
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3, 



here the fingers are allowed 



to bend along the gradient of the Laplacian field. The method of iterated conformal maps 
was applied to solve the growth problem in the geometry of the channel with two reflecting 
walls. We derived the Loewner equation for such a geometry and analyzed its analytical and 
numerical solutions in few simple cases. It was shown that the finger growth in the channel 
is qualitatively different from that in the unbounded space. In particular, the competition 
between the fingers is much stronger in the case of the channel and, for any positive value 
of the exponent 77, the only stable asymptotic situation is that of a single finger, which had 
outgrown the others and continues to grow with a constant velocity along the centerline of 
the channel, whereas the other fingers stop growing completely. 
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APPENDIX A: EXPRESSION FOR FINGER VELOCITY IN TERMS OF THE 

MAP ft 

In this appendix a sketch of the proof of Eq. (fTOjl is presented. This equation relates the 
finger velocity to //'(a) - the second derivative of ft mapping calculated at the image of the 
tip. To prove ( ITOll . it is convenient to introduce the analytic complex potential 6, such 
that u = ImB. Then, in the u plane, corresponding to the empty half-plane (or channel), 
the solution satisfying the boundary condition u{lj) = on the real axis is Q{uj) = u. 
The corresponding complex potential in the z plane is then simply Q{gt{z)) = gt{z). The 
derivative of the complex potential is directly connected to the gradient of u, in particular 

dgt 



dz 



\Vu\ 



(Al) 



Ignoring for the moment the singularity at the tip {z = 7), let us try to calculate the 
derivative of the g at that point 



dgt 

dz 



2=7 



lim 

<5^0 



6 



lim 



a + e{6) — a 



6 



e{6) 



lim 



(A2) 



where a + e is the image of the point 7 + 5 under gt. 
On the other hand 



6 



Ma + 6) - Ma) = f[{a)e + -r^ia)^ + ... 



(A3) 



The first term on the right hand side vanishes, since / has a local maximum at a, corre- 
sponding to the tip of the finger, thus e = (25//"(a))^/^ and the expression for the gradient 
takes form 



dgt 
dz 



z=7 



lim(5/;'(a)/2) 

0— >u 



-1/2 



(A4) 
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with the singularity at 5 = 0, as expected. Using the regularization described in Sec. HTl to 
remove the factor we get finally 



\n'ia)r/\ (A5) 



APPENDIX B: NUMERICAL METHOD 



In this appendix, we describe briefly the numerical method used to follow the evolution 



of the fingers. Instead of integrating the Loewner equation, as it is done e.g. in Ref. [47|, 
we obtain the mapping gt by direct iteration of the elementary slit mappings (which are 
of the form (fTTll for the half-plane and ( 1391 ) for the channel geometry). Since there are n 
fingers, each timestep involves the composition of n slit mappings, 0j, each characterized by 
a corresponding position of the pole, a^, and the growth factor di. In the case of the cylinder, 
the corresponding slit mapping may also be obtained from Eq (l39l) by a procedure described 
in Sec. |Vl for a finger with the pole at first the system of coordinates is translated by 
z ^ z — tti, then the slit mapping (l39l) is applied (the corresponding pole is now located at 
the origin, thus a = is to be put in Eq. f l38ll ) and the result is transformed back to the 
original coordinate system: z ^ z + ai. 

The calculations are somewhat complicated by the dependence of the growth factors di 
on the ft mapping, as given by Eq. ( l23l l. However, as before, this map may also be obtained 
by the composition of elementary mappings 0, which are the inverses of slit mappings (p, 
i.e. (j){(j){z)) = z. The second derivative of / may then be calculated either analytically, 
by differentiating the composition of all the slit mappings which make up ft, or by a direct 
numerical differentiation. The latter method is faster, whereas the former is more accurate 
which becomes important at the late stages of finger competition, where the growth factor 
of the losing finger becomes very small. 
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